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We introduce and simulate a growth model of the world-wide Web based on the dynamics of 
outgoing links that is motivated by the conduct of the agents in the real Web to update outgoing 
links (re)directing them towards constantly changing selected nodes. Emergent statistical correlation 
between the distributions of outgoing and incoming links is a key feature of the dynamics of the 
Web. The growth phase is characterized by temporal fractal structures which are manifested in the 
hierarchical organization of links. We obtain quantitative agreement with the recent empirical data 
in the real Web for the distributions of in- and out-links and for the size of connected component. 
In a fully grown network of N nodes we study the structure of connected clusters of nodes that are 
accessible along outgoing links from a randomly selected node. The distributions of size and depth 
of the connected clusters with a giant component exhibit supercritical behavior. By decreasing 
the control parameter — average fraction of updated and added links per time step — towards 
P C (N) < 10% the Web can resume a critical structure with no giant component in it. We find a 
different universality class when the updates of links are not allowed, i.e., for j3 = 0, corresponding 
to the network of science citations. 



I. INTRODUCTION 

Understanding the evolution of complex networks 
|l] [b| remains a challenging problem of modern statisti- 
cal physics. Modeling the dynamic structure of these net- 
works has significant practical applications, for instance 
in the design of search and control algorithms and in 
predicting emergence of new phenomena in information 
networks |^| . An evolving network consists of increasing 
number of nodes, which are connected by links accord- 
ing to certain rules specific to each network. Although 
the spatial distribution of nodes in a network is random, 
a different structure can be observed when the distribu- 
tion of links is considered. In this respect two major 
groups |)| are: (1) random networks, with the number 
of links attached to a node fluctuating around an aver- 
age value, and (2) networks with hierarchical structure 
of links, manifested in the power-law decay of the dis- 
tributions of node ranks. The structure of connections 
has immediate impact on the accessibility of nodes and 
on functional stability of the network. For instance, the 
hierarchically connected networks are shown ]To| to be 
robust to large failure rates, however, they are vulnera- 
ble to removal of a few 'key nodes'. The inhomogeneous 
structure of connections has been found in metabolic cy- 
cles pi, in various social networks 0], including the sci- 
ence citation network ||, and the Internet [Q and the 
world-wide Web §-§. 

The world-wide Web is not a separate physical net- 
works. Rather it is a subset of the Internet that uses 
a specific protocol (hypertext transfer protocol) and di- 
rected links to access a variety of data, representing 
an example of distributed client/server computing [ EJ . 
Based on these properties a hierarchical structure of links 
emerges that is peculiar to the world-wide Web. Recently 



large-scale measurements in the world-wide Web that in- 
volved 10 s nodes reveal || that: (a) Both the distribu- 
tions of incoming links and distribution of outgoing links 
follow a power-law behavior with different exponents; 
(b) The distribution of size of connected component — 
number of linked nodes in the web crawls — also exhibits 
scaling behavior; and (c) The Web shows a very intricate 
structure of connections, with a large strongly connected 
component in the center. 

Recently a great deal of effort was devoted to under- 
stand emergence of the scale- free structure of links in var- 
ious complex networks [p~p| p~3|] . An important ingredient 
of the dynamics — preferential attachment |Io| l was shown 
to lead to the power-law distribution for incoming links 
(the outgoing links are fixed by the rules of the model 
) . The model of preferential attachments for incoming 
links has been generalized to include initial attractiveness 
of nodes in Ref. pH , where the authors proved analyti- 
cally the existence of the scaling region for large evolution 
times. Recently it was demonstrated in Ref. 13 how a 



giant connected component can occur in the model of 
random directed graphs assuming arbitrary and statisti- 
cally independent distributions of incoming and outgoing 
links. The predictions of the random graph theory, how- 
ever, show substantial quantitative differences compared 
to the properties of the real world-wide Web Jll| . 

In this work we simulate growth and response of a di- 
rected graph with the dynamic rules that are motivated 
by prominent features of the world-wide Web. These dy- 
namic rules yield the statistically correlated distributions 
of outgoing and incoming links. Another important fea- 
ture of the model is that the links between pairs of nodes 
are not fixed in time, but may vary on the time scale of 
the network's evolution (updates of links) . This property, 
that may be found in some metabolic cycles, is not shared 
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with many other social networks, whose physical links are 
either fixed or vary on much slower time scale. We ar- 
gue that the (at present high) frequency of updates of 
the outgoing links, which is peculiar to the agents in the 
Web, is essential for the observed inhomogeneous struc- 
ture of connections. We find a consistent agreement with 
the recent empirical data on the real world-wide Web 
both in the structure of links and in the response of the 
network. 

The paper is organized as follows: In Sec. II we in- 
troduce the growth rules of the network and determine 
the emergent rank distributions of outgoing and incoming 
links. In Sec. Ill we study the temporal fractal structures 
in the growth phase of the network and determine the 
corresponding distributions of the first-return links. Sec. 
IV is devoted to investigation of the connected compo- 
nents on the fully grown network and fractal properties of 
noise. A short summary and the discussion of the results 
is given in Sec. V. 



II. GROWTH MODEL AND RANK 
DISTRIBUTIONS 

We suggest the following simplified set of rules to in- 
corporate basic features pertinent to the real world-wide 
Web: (i) Directed nature of linking. The world- 
wide Web represents a directed graph with nodes corre- 
sponding to Web pages, and arcs corresponding to hy- 
perlinks between pages |]|n|- (j) Growth and rear- 
rangements at unique time scale. At each time unit 
t a new node i = t is added to the network (growth) 
and a number M (t) of new links are distributed among 
the nodes (update of links) following two rules specified 
below. A fraction fo(t) = aM(t) of new links are out- 
going links from the new added node i = t, whereas the 
rest fi(t) = (1 — a)M(t) are the updated links at other 
nodes in the network. Hence, the relevant parameter in 
the model is the ratio of updated and added links at each 
time step, i.e., /3 = fi(t)/fo(t) = (1 — a)/a, which is in- 
dependent of the actual number of links M(t). Thus, in 
the model M(t) represents a net increment of the num- 
ber of outgoing links at time step t. Variations in M(t) 
are, in principle, not restricted, being caused by adding 
new connections between old nodes (resulting in the in- 
crease of the number of links), and/or by removing some 
earlier links (causing decrease of the number of links). 
We assume that variations in M(t) are such that we can 
define an average value M = M(t), which can be consid- 
ered as a constant in first approximation. In practice, the 
number of nodes and the number of links in the network 
increases with time, so that reasonable values for the av- 
erage M are positive. For consistency, we keep M = 1 
throughout this work. Different values of M do not af- 
fect the universal properties in the scaling region (i.e., 
for large evolution times). Taking much larger M val- 



ues, however, considerably increases computation time, 
(k) Preferential update (kl) and preferential at- 
tachment (k2). These properties represent a paradigm 
of social behavior, for instance preferential attachments 
are driven by "popularity" of a node, as discussed in ear- 
lier models of evolving social networks. In addition, here 
we have that not all nodes are getting updated at ev- 
ery time step, rather only a few nodes update at a time. 
Moreover, some of the nodes update outgoing links more 
frequently than others. This features can be formulated 
in terms of probabilities as follows: Where do the links 
come from? We assume that apart from the new node, 
updating at time t occurs with larger probability at most 
active nodes (preferential activity), i.e., an outgoing link 
from the node k < i appears with the probability 



Probl 



aM + q out (k,t) 
(1 + a)M * i 



(1) 



A new link goes to the node n (whose age is t — n) with 
the probability 



Prob2 = 



aM + q m (n, t) 



(1 + a)M * i ' 
a "popular" node attracts even more links. 



(2) 



i.e., a "popular" node attracts even more links. Here 
Qout(n,t) and qi n (n,t) are the dynamically varying num- 
ber of outgoing and incoming links, respectively, at the 
node n at current time t. It is assumed that at the time of 
addition of a node i to the network q ou t(i, i) — Qin(h *) = 
0. Note that the second rule in Eq. (||) is formally equiva- 
lent to the model of preferential attachment for incoming 



links |11|. However, in our model the role of the param- 



eter a in this equation is precisely determined through 
the rule in Eq. (|l|), which regulates outgoing links: (1) 
For a = 1/(1 + (3) < 1 a fraction (1 — a)M of outgoing 
links in Eq. ([j]) refers to updates at earlier nodes; (2) 
For a = 1 (equivalent to (3 = 0) no updates are allowed. 
New links originate only from the new added node; (3) 
Formally one can apply Eqs. (0^) for a strictly larger 
than unity, however, in this region (3 becomes negative, 
suggesting that the number of links in the network is not 
conserved. Therefore only the incoming links (^) can be 
counted correctly when a is in the interval oo > a > 1; 
(4) In the mathematical limit a — > oo, or in practical 
cases when a ^> q/M, the effects of the dynamical quan- 
tities q ou t(k,t) and qi n (n,t) become negligible, and we 
recover the case of fully random network. 

In what follows we would like to demonstrate that the 
features (i), (j), (kl) and (k2) explained above are essen- 
tial for the behavior observed in the Web. Earlier models 
fn^fnj concentrated only on the preferential attachment 
rule (k2), while largely neglecting properties (i), (j) and 
(kl). The most complete account of the preferential at- 
tachment model is given in Ref. |o]. On the other side, 
the random graph model [|l3| considered the directed na- 
ture of network (i), while neglecting the details of the 
linking rules (j), (kl) and (k2). 
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Using the rules in Eqs. we grow a large network 

of N = 10 6 nodes. We measure the ranking of nodes in 
this network by the distributions of outgoing links and 
incoming links, P(q ou t) and P(qi n ), respectively, which 
are shown in Fig. 1 for parameter (3 = 3, corresponding 
to a = 0.25 in Eqs. For comparison we also sim- 

ulate the distributions of outgoing and incoming links in 
the case of fully random directed network. For a range 
of values of their arguments the cumulative distributions 
for finite f3 values exhibit a power-law behavior according 
to 



(3) 



By increasing the parameter a in Eqs. (Q)-(g) in the range 
(0, 1), corresponding to decrease of the relative fraction 
of updated links (3 in the interval (oo,0), the slopes of 
the distributions increase. The scaling exponents n n 
and T out given in the inset to Fig. 1 are parametrized 
by a. Our numerical results for Tj n agree with the ex- 
act solution Jll]] Ti n = 2 + a for the incoming links. For 
the distribution of outgoing links no analytical results 
are available. Our numerical results in the inset to Fig. 
1 can be well approximated with the linear dependence 
r out ~ 2 + 3a. Therefore, the difference between the two 
exponents A = r out — Ti n « 2a persists for any finite 
< a < 1. Comparing these results with the available 
empirical data we estimate the parameter a = 0.22 ±0.1. 
Here we used the following values as the experimental 
estimates for the exponents Ti n — 2.16 ± 0.1, obtained 
as the average value from the data in [|| and last refer- 
ence in |10| , and T out = 2.62 ± 0.1 from data in ||, by 
fitting only the straight part of the curves while avoid- 
ing the noisy tails and the curvature at small q. It is 
clear that these values are not definitive. More careful 
fits of the data and further measurements are necessary 
in order to reduce the experimental error bars. Hence, 
the approximative value of the control parameter /3 in 
the current state of the Web is (3 w 3, i.e., in the average 
three updated links come to one added link at each evo- 
lution step. We now discuss the properties of the network 
for this particular value of the control parameter (3. 



III. DYNAMIC FRACTAL STRUCTURES 

In the dynamical systems with large-scale organization 
the occurrence of algebraically decaying distributions can 
be linked to the formation of certain fractal structures 
M| in the system. In the case of complex evolving net- 
works, however, the existence of such structures is less 
clear, due to the spatial randomness of the graph. To 
show that the power-law decay of the distributions of 
node ranks, shown in Fig. 1, is associated with certain 
fractal structures on the network we examine the tem- 
poral pattern of linking between nodes. A part of the 



pattern is shown in Fig. 2 (top panel). Clustering in the 
patterns indicate frequently active nodes, whereas voids 
of various sizes indicate lesser activity at corresponding 
nodes. Visually heterogeneous picture suggests that both 
pattern of origins of links and pattern of targets have 
fractal structure. The fractality of these patterns can be 
quantified, for instance, by computing the distribution of 
time intervals between two successive linking at a given 
node (return time). The distribution of time intervals At 
between two successive linking to a selected node in the 
network P(At) measured for the network of 10, 000 nodes 
and the parameter (3 = 3 is given in Fig. 2 (lower panel). 
The algebraic decay of the distribution is a signature of 
the underlying fractal structure of the pattern. The dis- 
tribution of time intervals between two successive links 
from a selected node also shows a power-law tail for large 
At. At small intervals At the pattern is nearly random. 
The two distributions coincide, indicating mutual corre- 
lations in the patterns, for time intervals At > 100 (that 
are accessible for large evolution times t). To emphasize 
the relation to Fig. 1 we recall that the number of links of 
either kind at a node accumulates with time. Hence, the 
regions of large time intervals between successive linking 
at a node, At in Fig. 2, correspond to small number of 
links at that node, q in Fig. 1, and vice versa. 



IV. DYNAMIC CRITICAL STATES 

The scaling behavior of the rank distributions can be 
examined in terms of the potential dynamic critical states 
by measuring the statistics of triggered avalanches (l^] on 
the network. We consider an analog of the avalanche size 
in the networks: the size of a cluster of nodes which are 
physically accessible along directed links starting from a 
random node in the network. A network is grown using 
the rules in Eqs. (@)-(|). We then select a random node 
and make a list of nodes which are connected by outgoing 
links from that node. In the next step we make a new list 
following outgoing links from the nodes on the previous 
list, keeping only the nodes which are different in the two 
lists, and so on. The process ends up when no new nodes 
can be reached — the list is empty. The number of differ- 
ent lists before an empty list occurs can be recognized as 
the depth of the connected component. Technically, we 
preserve ranks qi n and q ou t for each node from the growth 
phase and, in the case of small networks, the exact physi- 
cal links between the nodes. In large networks only ranks 
are preserved and links are searched using the rule in Eq. 
(|l|). For a large ensemble of networks the results are ex- 
pected to be statistically the same. In this way the size 
of a cluster of linked nodes represents a response of the 
system on the random excitation. The process is remi- 
niscent of the invasion percolation on a random graph, 
where a new edge j, which already does not belong to 
the graph, is invaded along an outgoing hyperlink i — » j 
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from the node i if the node rank q ou t(i) exceeds unity. 
The probability to add the edge j to the graph is given 
by the preference rule in Eq. (|]) . In a more familiar case 
of the invasion percolation in physical wetting on the lat- 
tice [jlTj that probability is given by the least resistance 
rule. To our knowledge the problem of invasion percola- 
tion on random graphs has not been studied so far pc( |. 

The differential distribution of size of connected clus- 
ters of nodes is given in Fig. 3 for different values of the 
parameter (3. The distribution of depth of the same clus- 
ters is given in the inset to Fig. 3. In the case (3 = 3 we 
find that small clusters follow a power-law distribution 
with the exponent t s = 2.79 ± 0.12, comparable with the 
measured [pi value t s — 2.52. (Note that in the case of 
Internet the corresponding exponent is estimated (J] as 
Tg = 1.90.) In Fig. 3 the occurrence of the peak in the 
distributions at large clusters So « 10 2 indicates that a 
large number of the avalanches have approximately same 
size. This is a signature of the existence of a large subset 
of strongly connected nodes (a giant component): Once 
a 'crawl' enters the subset of strongly connected nodes, it 
explores it entirely. By taking a larger number of nodes 
N the position of the peak moves towards larger values 
approximately as So ~ N/7. 

Before discussing the distributions in Fig. 3, we demon- 
strate that the network with the hierarchical structure of 
links is characterized by a fractal noise. In the search of 
a connected cluster described above we examine detailed 
variation of the number of nodes added to the cluster 
at each step of investigation. We obtain the noisy sig- 
nal shown Fig. 4. Properties of the signal vary with the 
control parameter of the dynamics /3. Considering the 
number of steps as a total elapsed time of investigation, 
the Fourier spectrum of the signal is shown in the top 
panel in Fig. 4. It shows a region of correlated behav- 
ior with power-law decay between the upper cut off at 
high frequencies and lower cut off (due to finite size of 
the network). In addition, the spectrum in the case of 
supercritical network exhibits a peak at a characteristic 
frequency /o, which is absent in the case of the critical 
parameter (3 = (3c (N) (cf. Fig. 4). 



V. DISCUSSION AND CONCLUSIONS 



In the theory of critical states |1^J19|], appearance of 
the peak in the distribution of avalanches, such as the 
distributions shown in Fig. 3, indicates that the system 
is supercritical. In the critical state the absence of any 
characteristic scale is manifested in the purely algebraic 
decay of the distribution of cluster size until a cut-off, 
which depends on the size of the network N. On the 
other hand, the form of distribution in Fig. 3 suggests 
that a critical point exists that can be reached by vary- 
ing a relevant parameter of the dynamics |ll| . In fact, by 
decreasing the parameter (3 we see in Fig. 3 that the slope 



t s decreases and the peak eventually disappears at a crit- 
ical value |Q of the parameter /3 = /3 C ~ 0.081. Beyond 
the critical value of the control parameter the scaling be- 
havior of the size and depth of connected components is 
entirely lost. In the critical state the scaling exponents 
are (cf. inset to Fig. 1): Ti n k, 2.925, t s = Ta = 1. 
The distribution of outgoing links shows very sharp de- 
cay, the exponent r out is difficult to measure. It should 
be stressed that the absence of the peak in the critical 
state suggests that no giant component can be formed. 
Rather, the network consists here of many small groups 
of well interlinked nodes and percolating directed links 
between these groups. 

It should be stressed that the special limit of our model 
(3 = (i.e., a = 1), corresponding to the original model 
of preferential attraction of incoming links proposed by 
Barabasi, Albert and Jeong fl0|| , belongs to an entirely 
different class of scaling behavior as regards the distribu- 
tion of size and depth of connected clusters: these dis- 
tributions are flat between lower and upper cutoffs (see 
Fig. 3). The only distribution which enjoys a power law 
in this limit is the distribution of incoming links, which 
has the exact |11 exponent Tj„ = 3 (see also inset to Fig. 
1). In this case ([3 — 0) updates at already existing nodes 
are no longer possible. This feature — freezing of the out- 
going links — as well as quantitative disagreement in the 
exponent 7 for distribution of incoming links, makes the 
model of Ref. |L0| inappropriate as a model of the world- 
wide Web dynamics. A realistic network with frozen out- 
going links is the network of scientific citations Q , where 
physical links, corresponding to the cited references in 
already published papers, remain fixed in time. 

In conclusion, we have demonstrated that certain 
salient features of the dynamics of the world-wide Web 
require more careful modeling, compared to models of 
a generic complex evolving network with preferential at- 
traction of links |iof and models of random graphs with 
rewiring fl3|| . The dynamic structure and functioning 
of the world-wide Web is deeply rooted in the activity 
of the agents who are creating the outgoing links. The 
hierarchical structure of outgoing links, which is docu- 
mented by measurements in the real world-wide Web [^) , 
is related to the bias activity of the agents in the rule 
(|l|) of our model. A random selection of the active agent 
(see curve (c) in Fig. 1) fails to describe the distribution 
of outgoing links in the real Web. Temporal variation 
of the outgoing links inside the network, i.e., updates of 
links, has twofold consequences on the global structure of 
the network: (i) When the updates of the outgoing links 
are allowed (i.e., when f3 strictly larger than zero), the 
structure of incoming links qualitatively changes, the ex- 
ponent Ti n (P) < 3. The distribution of incoming links in 
our model coincides with the analytical results of Doro- 
govtsev et al. |ll]] in the scaling region. In addition, the 
distribution of size and depth of connected clusters be- 
comes hierarchical when (3 > 0, in agreement with results 
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found in the real Web, which is not the case when no up- 
dates are allowed (/3 = 0). (ii) Varying the frequency 
of updates (3 implies changes in the Web structure, both 
in the outgoing and incoming links. This has immediate 
impact on the accessibility of nodes. The correlations 
between the outgoing and incoming links suggests that 
the local structure of the network is qualitatively different 
compared to the case without updates. Here we did not 
study in detail the probability that two nodes are linked 
(clustering coefficient). Rather we studied the physical 
properties of the network which make the background 
of the observed behavior: fractal temporal patterns and 
number of nodes that join the cluster at one time unit. 

By comparison of the simulated results and the data 
obtained in the real world-wide Web we obtained a sys- 
tematic agreement both for the distributions of outgoing 
and incoming links and for the size of connected com- 
ponents when a single control parameter is fixed to 
= 3 within error bars of the data. To our knowledge 
no previous model claimed to describe the world-wide 
Web achieved such degree of consistency. This makes us 
believe that the present simplified model takes into ac- 
count some basic features of the dynamics of Web. It 
would be interesting to estimate (3 by directly measuring 
the average number of updated links in the world-wide 
Web relative to the number of added links originating 
from each added node. 

The structure of links shown by the distributions in 
Figs. 1 and 2 is closely related to the evolution of the 
number of links at a given node and to other local prop- 
erties of the network. A detailed study of these prop- 
erties within the present model requires additional work 
that remains to be done in the future. Here we expect, 
based on the analytical results for the incoming links of 
Rcf . jll| , that in our model the average number of links 

< qin{i,t) > at a node i will decay in time t > i as 

< q(i,t) >~ (i/i)~ 7 "*, where the exponent ji n is given 
by the exact scaling relation [|ll| j in — l/(r in — 1), lead- 
ing to ji n = 1/(1 + a) ss 0.86. Assuming that the density 
of the outgoing links at a given node also exhibits scaling 
behavior, as Figs. 1-2 suggest, we can predict a slower de- 
cay for the average number of outgoing links at a given 
node. The expected exponent is j ou t ~ 1/(1 + 3a) ~ 0.6. 

The simulation results that we reported here suggest 
that the emergent structure of links in the world-wide 
Web is strongly related to the updating policy of the 
agents: who updates and how often. It remains to under- 
stand the potentially more intricate reverse effect: how 
the amount of information currently stored in the Web 
influences the conduct of the agents, with implications 
on self-tuning of the control parameter (3. In our model 
the structure of the network at j3 > f3 c (N) appears to 
be supercritical, where (3 C {N) < 0.1 in a large network. 
Therefore, it is plausible to expect that the evolutionary 
selected values of (3 will be much smaller than the cur- 
rently estimated value (3 « 3, if the scenario of self-tuning 



of the control parameter is active in the real Web. In this 
respect the current state of the real Web can be regarded 
as a transient rather than a stationary state. 

Finally, our results support the conclusion that in the 
directed graphs two growth rules are necessary to de- 
scribe the dynamics of the outgoing and incoming links, 
respectively. In the case of the world-wide Web the statis- 
tically correlated distributions of outgoing and incoming 
links appear as a fundamental feature of the evolution of 
the Web. This conclusion may serve as a starting point 
for the future modeling of the real world-wide Web in 
terms of the master equations. 
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FIG. 2. Top panel: Part of the temporal pattern of linking 
in the growth phase of network with TV = 10, 000 nodes. One 
link per time step is considered. Dots represent nodes from 
which the link originates, whereas crosses are target nodes. 
Lower panel: Distribution of time intervals At between two 
successive linking (a) from a given node, and (b) to a given 
node. The two patterns become correlated for the time in- 
tervals At in the range 100 < At < 3, 000, corresponding to 
the scaling region of the rank distributions 2000 > q > 200 
in Fig. 1. Distributions are averaged over 100 samples and 
logarithmically binned. Both distributions are normalized to 
the total number of time steps. Note that the events with 
At = that correspond to the outgoing links from new added 
nodes, contributing to the distribution (a) , are not shown. 



FIG. 1. Cumulative distributions of outgoing links (a) and 
incoming links (b) for the network of TV = 1,000,000 nodes 
and average ratio /3 = 3 of updated respective to added links 
per time unit. For comparison we have included the corre- 
sponding distributions (c) and (d) in the case of fully random 
directed graph. Fitted slopes of the straight sections of the 
curves (a) and (b) are compatible with the scaling exponents 
in Eq. (3) as r ou t — 2.75 and Ti„ = 2.25, respectively. Inset: 
Variation of the scaling exponents T out — 1 (diamonds) and 
Tin — 1 (triangles) with the control parameter a = l/(/9 + 1) 
in the physical range (0, 1). Dotted line: exact solution for 
the case of incoming links from Ref. [11]. 
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FIG. 3. Differential distribution of size of connected clus- 
ters in the network grown with the rules of Eqs. (|l|)-(^|) for 
different values of the control parameter [3 = 3 (diamonds), 
1 (circles), 0.081 (triangles), and zero (crosses). We employ 
1, 000 avalanches in each of 100 samples of the network with 
N = 1,000 nodes. Inset: Corresponding distributions of 
depth of the connected clusters. 
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FIG. 4. (Lower panel) Number of nodes n(t) added to 
a connected cluster at one investigation step vs cumulative 
number of steps t. Size of one connected cluster is repre- 
sented by the surface enclosed between two consecutive drops 
of the signal n(t) to the base line n(t) = 1. Parameters 
are: N — 1,000 and j3 = 1. (Top panel) Fourier spec- 
trum of the same signal. Also shown is the spectrum of the 
corresponding signal at the critical value of the parameter 
f3 c {N) = 0.081. Data are logarithmically binned. Straight 
lines are power-law fits with the slopes <j> — 1.45 ± 0.06 (for 
13 = 1) and </> = 1.01 ± 0.06 (for = 0.081). 



7 



